jj — H Lamb shift: dispersing the nucleon-excitation uncertainty 
with a finite energy sum rule. 
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We assess the two-photon exchange contribution to the Lamb shift in muonic hydrogen with 
forward dispersion relations. The subtraction constant T(0,Q 2 ) that is necessary for a dispersive 
evaluation of the forward doubly-virtual Compton amplitude, through a finite energy sum rule, is 
related to the fixed J = pole generalized to the case of virtual photons. We evaluated this sum 
rule using excellent virtual photoabsorption data that are available. We find that the "proton po- 
larizability correction" to the Lamb shift in muonic hydrogen is —(40 ± 5)fieV. We conclude that 
nucleon structure-dependent uncertainty by itself is unlikely to resolve the large (300/ieV) discrep- 
ancy between direct measurement of the Lamb shift in \iH and expectations based on conventional 
Hydrogen measurements. 



I. INTRODUCTION 

An ongoing controversy surrounding the proton size 
originates from the large discrepancy between the re- 
cent measurement of the Lamb shift in muonic hydrogen 
and earlier measurement based on conventional hydro- 
gen as well as electron scattering (see for example the 
review [T|). The advantage of using the muonic hydro- 
gen over the conventional is that due to a larger reduced 
mass the Lamb shift in the former is by an order of mag- 
nitude more sensitive to the proton radius. The Lamb 
shift AE2P-2S m muonic hydrogen depends on the pro- 
ton charge radius, Re through [2] 



AE 2P -2S = 205.93(1) meV - ~aen<% s (0)R% + 0(a 



(1) 



where the wave function at origin is given by 
(am r ) 3 /8ir, a = e 2 /Air is the fine structure constant, and 
m r = m p M/(m M + M) is the reduced mass with m M , M 
the muon and proton masses, respectively. The value of 
the Lamb shift predicted using Re quoted by the Com- 
mittee on Data for Science and Technology (CODATA) 

m 



R E = 0.8768 (69)fm, 



(2) 



that is based primarily on the electronic Hydrogen Lamb 
shift measurement, or on the value extracted from the 
most recent electron scattering data [3], 



R E = 0.879 (8) fm, 



(3) 



differs by 5er from the measurement of the muonic hydro- 
gen Lamb shift by Pohl et al. [S]. The later requires a 
significantly smaller charge radius, 



R E = 0.84184 (67) fm. 
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(4) 



In terms of the Lamb shift, the discrepancy amounts to 
some 300/xeV that by far exceeds the experimental sen- 
sitivity of the muonic experiment [5 J. The first term in 
Eq. (fTl), that represents, up to 0(a 5 ), all QED effects 
associated with the leptonic current is almost three or- 
ders of magnitude larger than the observed discrepancy. 
This may lead to a conclusion that a slight adjustment in 
one of those terms could resolve the whole puzzle. These 
higher-order QED corrections, however, have been known 
for a long time and are well established. The reader 
is referred to three recent reviews which assess the full 
body of the relevant QED corrections, [2l|6l[7j)- A non- 
perturbative numeric evaluation is also available [5] and 
yields a similar result, and so does the analysis based 
in effective non-relelativistic expansion of QED fTU] , 
An exotic possibility is a substantial non-universality of 
lepton-proton interaction, which has not been observed 
before but a more plausible explanation is that higher 
order terms in the expansion in a is responsible for the 
discrepancy. Since QED corrections have a solid found- 
ing, attention has been focused on higher-order, nucleon 
structure-dependent effects. To lowest order, 0(a 5 ) these 
arise through a two-photon exchange process and poten- 
tially bear significant uncertainty because they involve 
the complete nucleon excitation spectrum. 



In Section [TTJ we assess this two-photon exchange con- 
tribution to the Lamb shift using forward dispersion re- 
lations. Section (TTTJ deals with the novel feature of our 
approach, were we use the finite energy sum rule (FESR) 
to relate the value of the subtraction function that arises 
in the dispersive calculation to the contribution from the 



fixed J = Rcgge pole. Section IV is dedicated to the 
numerical analysis. Discussion of the results and com- 
parison with the existing calculations is summarized in 
Section [V] 
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II. DISPERSION RELATIONS FOR COMPTON 
SCATTERING 

The 0(a 5 ) contribution to Lamb shift sensitive to pro- 
ton structure enters through the matrix element of the 
two-photon exchange (TPE) between the lepton and nu- 
cleon integrated over the atomic wave function. This 
can be seen as the virtual excitation and de-excitation 
of the proton by the successive photons, and thus all the 
complexity of the excited nucleon states is affecting a 
precision atomic physics computation. Taking the stan- 
dard approach for computing bound state corrections in 
atomic physics which express nucleon current effects in 
terms of the atomic wave function at the origin the TPE 
contribution to the Lamb shift is then given by [TTJ [12] 

F = A M°) 2f A {q 2 + 2 V 2 )T 1 -{q 2 - V 2 )T 2 
2m; J (2tt)4 q%q 2 /2 mi ) 2 - u 2 ] ' 

(5) 

where m,, i = e, \x is the lepton mass in conventional 
and muonic hydrogen, respectively. The scalar functions 
Ti.2 = Txpiy, q 2 ) with v — (pq)/M, are the standard am- 
plitudes that parametrize the spin-independent hadronic 
tensor for doubly virtual forward Compton scattering 
7* (q) + N(p) — > 7*(g) + N(p), and are given by 

= -i— / d 4 xe lqx (N\T[J^{x) i r(0)}\N) 
ottM J 

= (-.^ ,y + ^r-) T i(^ 2 ) ( 6 ) 



relation [15] , 

oo 

oo 

ReT 2 (,,Q 2 ) = fvf (8) 
2n J (v' z — v l ) 

While this suffices to reconstruct T 2 from knowledge of 
the dispersive part, T\ requires an additional input in the 
form of a subtraction constant at each Q 2 , i.e the func- 
tion 2i(0, Q 2 ). This is due to divergence of the unsub- 
tracted dispersive integral at large energies as dictated by 
the high energy asymptotic properties of the F\ structure 
function. At the real photon point Q 2 = 0, the subtrac- 
tion term is fixed by the well-known Thomson-scattering 
limit, 7i(0, 0) = —a/M. For virtual photons however, 
existing estimates carry large uncertainties. They are 
based on the not so well determined polarizability and 
the Q 2 dependence of elastic form factors. 

The Fi structure functions measured with virtual pho- 
tons receive a contribution from the single nucleon pole 
(Born terms) at v tr = vn = ±Q 2 /2M , and from the uni- 
tarity cut due to opening of particle production thresh- 
olds which start with pion production at v tr — v^{Q 2 ) — 
±[(M + m^) 2 - M 2 + Q 2 ]/2M (with being the pion 
mass). Following [TT], we divide the contribution to the 
Lamb shift into three physically distinct terms that orig- 
inate from the subtraction term Xi(0, Q 2 ), the nucleon 
pole and finally all excited intermediate states that may 
couple to 7 AT, respectively 

AE = AE subt + AE el + AE mel . (9) 



The hadronic tensor can be measured in a restricted 
kinematic range of the variables v and Q 2 and needs to 
be extrapolated outside the physical range to compute 
the integral in Eq. (J5j|. The extrapolation is based on 
analytical continuation. Specifically, the functions Ti )2 
are discontinuous along the real axis in the complex en- 
ergy plane v with the discontinuity, which is equal to the 
imaginary part, related to the inclusive cross section 



with 



ImTi(i/,<2 2 ) 
lvoT 2 {v, q 2 ) 



AM 



(7) 



As customary in dispersive approaches, we make use of 
the complex v — (s — u)/{AM) plane. Since this vari- 
able is crossing-symmetric, upon applying Cauchy's the- 
orem, the left and right cut can be combined in the same 
integral, yielding a relatively simple forward dispersion 



AE' 



AE el = 



Ti(0,Q 2 ) 




miM 



Mv 

- 7 y^l2{T,T l )F 2 {v, 



(10) 



n 



V(4mf 



7(4Af 2 ), r = v 2 /Q 2 , and the 
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auxiliary functions denned by 



71 (r) = (1-2t) v / 1 + t + 2t 3/2 

72 (r) EE (l + T )3/2_ r 3/2_3^ 



7i (t, 7 "/) 



^71 (~n) - V^7i(t) 

Tl - T 



~ fV - ^ - 1 f 72 (t) 72 (tQ 



(11) 



Note that generally, besides the integral over the union 
continuum that is represented in the above equations, a 
sum over the discrete spectrum must be taken. The lat- 
ter contributes to the Lamb shift at order 0(a 6 ) and is 
dropped from our considerations. Using these formulae, 
in [TT] the inelastic contribution, AE mel was evaluated 
using the photo-absorption cross section parametriza- 
tion of |14j for the resonance region complemented with 
the high energy paramctrization of 15J. Their elastic 
(nucleon-pole) contribution, AE el was computed using 
three different phenomenological parametrizations of nu- 
cleon electromagnetic form factors [H [THl Q2] . Here we 
also give an independent evaluation of the two contri- 
butions. For AE mel we use a recent parametrization 
of inclusive structure functions [T5] that also uses the 
parametrization of the resonance region from |14) but it 
uses a modified Regge-inspired background that is fitted 
to the total photoabsorption cross section of 19J. The 
Q 2 -dependence is introduced as in [20]. For E el , we use 
the parametrization from |16) to finally obtain 



AE e 



-30.1 neV, AE V 



-13.0 £teV (12) 



Within errors these agree with computation in the origi- 
nal analysis of [TT] 



AE el = -29.5 ± 1.3 fieV, AE mel = -12.7 ± 0.5 fieV 



(13) 



III. EVALUATION OF THE SUBTRACTION 
TERM 

A. Finite energy sum rules 

While previous analyses concentrate on the low energy 
constraints for the subtraction term, here we focus on 
implications of the high energy behavior for constrain- 
ing the subtractions. This is done by exploiting the fi- 
nite energy sum rule (FESR) for the Compton amplitude. 
The subtraction term in the dispersion relation (DR) for 
Ti arises because the high-energy photo absorption cross 
section does not vanish asymptotically. It can be well 
described by a Regge-theory inspired parametrization 



a p i 



Cfl(0) 



a R -l 



(14) 



with the effective Pomeron and leading Regge trajectory 
intercepts given by ap = 1.097 and ap = 0.5, respec- 
tively. The contributing to the cross section is deter- 
mined by cp(0) = 68.0 ±0.2^b and cr(0) = 99.0±1.2^b, 
with Vq = 1 GeV. 

The corresponding contribution to the Compton am- 
plitude T\ of this Regge part is given by 



(15) 



ReT^i/,0) = 



2tt 2 



dp 



,'2 



Following [22 , we write a dispersion relation for the 
difference, T x - T?, 



Re Ti (v, 0) - Re T^(v, 0) = 
M 2tt 2 



(16) 



With the large-f tail thus removed, the dispersion inte- 
gral on the right hand side of Eq. ( 16 ) is dominated by 
energies below a scale N = O(v ) which is discussed be- 
low. Removal of the asymptotic contribution from the 
dispersive integral introduces a new subtraction, Coo de- 
fined by, 



Coo(0) = [ReT x (i/,0) - ReT/V,0)] 



(17) 



With the help of currently available high energy data, 
Cqo (0) has recently been determined with high accuracy 
[23"] and it follows from Eq. ( 16 ) that it is related to the 
high energy parameters by 



Oo(0) 



+ 



a 
2tt 2 



E 

,=P,R 



1 

2^2 

c.-(0) 
ou 



N 



dv'a T (v',Q) 



(18) 



The resonance contribution given by the integral over the 
photoabsorption cross section is well established and can 
be readily evaluated from the low energy data. The pa- 
rameter N defines the lowest photon energy above which 
Regge parametrization suffices to describe the data, 
which in the analysis of [23 was taken to be 2 GeV. From 
this analysis it follows that Coo(0) = (—0.72 ± 0.35) fib 
GeV. 

For our application to muonic hydrogen we need to 
generalize the above, real Compton amplitude dispersion 
relation to the virtual photon case. Using the relation 



Fxiv, Q 2 



Mv{\ - x) 



(19) 



where x = Q 2 / (2Mu), we may write 



T^Q 2 ) =T!(0,Q 2 



dv'F^v'.Q 2 ) 
■InM J v'{v' 2 -v 2 ) 



(20) 
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In analogy to the real photon case we introduce the 
Regge-theory motivated representation for the high- 
energy data valid for v > N(Q 2 ), 



p rpRf n 2, ^ [°° , Ff(l/', Q 2 ) 



with 



Mi/ 



i=P,R 



(22) 



The generalization of Eq. \22\ is not unique since 
in principle vq and on might be made Q 2 -dependent. 
These eventual Q 2 -dependences for low Q 2 < 1 GeV 2 
that are of interest here can however be absorbed in 
Ci{Q 2 ) without loss of generality. The coefficients Ci(Q 2 ) 
must reduce to those found for real photons at Q 2 = 
that are listed below Eq. (14). Their Q 2 dependence, 
and that of N(Q 2 ), is obtained by matching the Regge- 
parametrization of Eq. (22) and F\(v, Q 2 ) defined by 
Eq. (Tl9| For v > N(Q 2 ) and moderate Q 2 < 1 GeV 2 , 



we obtain 



cp(Q 2 ) = c P (0) 

cr(Q 2 ) = c*(0)-(20±10)/xb 



Q 

GeV 



and 



N(Q 2 ) 



>GeV 



2M ! 



(23) 



(24) 



Note that the presence of the factor 1—x = \—Q 2 /2M v 
in the relation between ut and F\, Eq. (19) requires a 



value of N(Q 2 ) larger than that found for real photons 
iV(0). In any case, the resulting FESR will not be sen- 
sitive to the value of N, as long as the Regge amplitude 
correctly represents the data for all v > N. The values 
cp(0), cr(0) are fixed by very precise fit to real photoab- 
sorption data, and cp(Q 2 ) is moreover fixed to its real 
photon value (for low Q 2 < 1 GeV 2 only) to ensure that 
asymptotically <tt — &t vanishes, the assumption that is 
crucial for the FESR method. This effectively leaves the 
Q 2 -slope of the coefficient cr(Q 2 ) taken as a linear func- 
tion the only parameter that has an uncertainty, and we 
assign a generous 50% uncertainty thereto. The analog 
of Eq. (fisb at finite Q 2 , 



C^Q 2 ) = [ReT 1 ( i y,Q 2 )-ReT 1 fl ( i /,Q 2 )]| 
satisfies now 



(25) 



N(Q ) 

- r dv 



C^Q 2 ) ^ l( 0,Q 2 )-— ; ; 

c t (Q 2 ) ( N(Q 2 ) 



TF^v'iQ 2 ) 



(26) 




FIG. 1: The residual term for the high energy Compton 
amplitude corresponding to a fixed pole at J = complex 
angular momenttum plane. It corresponds to Compton scat- 
tering on a pointlike quark at instant light-cone time. 




FIG. 2: Regge exchanges in the t-channel dominate the high- 
energy part of the Compton amplitude. 



It is expected, that at high Q 2 Coo(Q 2 ) is finite and rep- 
resents a light-cone instantaneous, two-photon interac- 
tion on a point-like quark |25j . as depicted in figur e [T[ 
This causes no problem in the first of equations ( 10 ) 
for E subt that is convergent upon substitution of a con- 
stant contribution to T^O, Q 2 )- The constant C OCJ (Q 2 ) 



is related to the virtual Compton amplitude Ti(0,Q 2 ) 
through Eq. ( 26 1 and enters the Lamb shift though E subt . 



To evaluate the integral on the right hand side of 
Eq. ( 26 1 we need a parametrization of the virtual photon- 



proton cross section to substitute in Eq. (19), for which 
we use the form obtained in [T5] fits to electroproduction 
data, 

<tt(W 2 ,Q 2 ) = ]TBIF a (FV 2 )F 2 (Q 2 ) (27) 

a? ot (W 2 ,0)F B (Q 2 ). 



In the first term the summation runs over nucleon reso- 
nances with BW standing for a Breit-Wigner propagator, 
BW a (W 2 ), and electromagnetic transition form factors 
given by F a (Q 2 ). The second term represents a smooth 
background. Expressing Ti(0, Q 2 ) in terms of the J = 



5 




FIG. 3: The low and intermediate energy region is described 
by a sum over a s-channel resonances that are photoexcita- 
tions of the nucleon. 



pole contribution, C OD (Q 2 ) yields, 



2tt 



N(Q 2 ) 



2Mv 



(28) 



MQ 2 ) 



which is the main result of this paper. It expresses the 
low-energy function Ti(0, Q 2 ) that enters the Lamb shift 
through E subt in Eq. (19]) in terms of three distinct con- 
tributions with clear physical interpretation, which are 
diagrammatically shown in figures [TJ [2] and |] The last 
two are the i-channel Regge exchanges and s-channel res- 
onance contributions; the split between the two is deter- 
mined by N(Q 2 ). The first term is the J = 0, fixed-pole 
contribution to virtual Compton scattering Cx>(Q 2 ) [24] 
to which we now turn our attention to. 



B. Analysis of the fixed pole 

The J = fixed pole in Compton scattering was in- 
troduced in [55] and studied in phenomenological models 
e.g. in \H\ [25] • Such an s and t independent 

contribution has been analyzed in the kinematic region 

-t ^> M$f and the existing 



where both —t, s are large, s, 
data in this region [29 , 30J supports existence of the fixed 
pole. 

For real Compton s catt ering Coo (0) was determined in 
[2"3"| . however, in Eq.(28) is evaluated at finite Q 2 . 
Theory suggests that at asymptotic Q 2 , C 00 (Q 2 ) is con- 
stant |25j . but this has not been experimentally estab- 
lished; it might be so in the future with the help of the 
Deeply Virtual Compton Scattering program at Jeffer- 
son lab. To allow for the possibility of a Q 2 dependence, 
we subtract Eq. ( 18 ) (real FESR) from Eq. 



FESR), and changing the integration variab 



28) (virtual 



e from v to 



oj = v- Q 2 /2M, obtain 



T 1 (0,Q 2 ) = -— + [C^Q 2 ) - ^(O)] 



(29) 



1 

~2^2 
~2tt 2 



N(0) 



dw 



MO) 

E 

=P,R 



Ul 



-<jt(w, Q 2 ) - o-t(w,0) 



2M 



Ci (0){N(0)\ ai c t (Q 2 ) fN(Q 2 )\ a > 



Fb(Q 2 ) 
(30) 



This is a rigorous representation of the subtraction term 
in the virtual Compton amplitude. If the fixed pole were 
Q 2 independent, as suggested by [25], Coo would drop 
out of this equation. Since this is not established experi- 
mentally, we also provide an order of magnitude estimate 
under the assumption that C OCJ (Q 2 ) falls with Q 2 . 

For the estimates of the uncertainty associated with 
Coo(Q 2 ) — Coo(0) we use a parametrization 



Coo(Q 2 )-Coo(0) 



A 2 + Q 2 



[CooM-Coo(O)] , (31) 



with a typical scale A = 1 GeV and Coo(oo) = 0. 



IV. NUMERICAL ANALYSIS 



If we substitute Eq. ( 30 1 in the expression for E in 



Eq. (10) we see that the result is IR divergent. This is 



due to the Thomson term, Ti(0, 0) = —-£?• Physically, 
it corresponds to exchange of soft Coulomb photons that 
is already taken into account at the level of atomic wave 
functions, and has to be subtracted in order to avoid 
double-counting. We are left with the following conver- 
gent integral to be evaluated: 



oo 

AE subt = 4a0 2 (O) JdQj^n) 



The contribution from 7i(0, Q 2 ) to the Lamb shift can 
be written as a sum of several terms, 



AE 



AEl es + AE Back + AE Re " e (33) 



We evaluated the respective integrals in Eq. ( 32 ) numer- 



ically. Below, we quote the individual contributions from 
each of the well-established resonances, the non-resonant 
background, and the Regge part, respectively, 



G 



A# A (i232) = (0.95 ± 0.09) /^eV 
A£ Sll (i535) = (-4.02 ± 3.14) M eV 
AE Dl3{1520) - (0.41 ± 0.09) /xeV 
(-0.23 ±0.16) neV 
(-0.32 ±0.06) /xeV 
(0.10 ± 0.02) ^eV 
(-0.76 ±0.26) /xeV 
A£ ,«ac« = (_29.34±2.93)^eV 

AE Regge = ( 36 . 5 5 ± 1 . 6 ) ^eV , 



AE, 
AE 
AE_ 
AE 



Sn(1665) 

Fib (1680) 

Pi i (1440) 

F 37 (1950) 
Back 



(34) 



Adding the above contributions to the subtraction term, 
AE subt = ^3 3 ± 4 6 ^ ^ eV ( 35 ) 

It can be noted that there are strong cancellations be- 
tween various terms. The size of the correction is al- 
most entirely given by the sum of three contributions, 
AE Reg 9 e AE Back and A£ 5ii(1535) . To discuss the un- 
certainty it thus suffices to constrain the uncertainty in 
these three contributions. Regge and background contri- 
butions are large, opposite in size and cancel to about 
80%. The background contribution is obtained from a 
fit to excellent experimental data over a wide range of 
W 2 ,Q 2 (see Ref.[T4] for a full list of references) and a 
relative uncertainty of 10% is reasonable. The Regge 
contribution is related to the background since they are 
constructed to coincide at high energies, and assigning 
an extra uncertainty here would lead to double count- 
ing. We assign a 50% uncertainty on the Q 2 -slope of the 
Reggeon strength cr(Q 2 ). For the resonances, we assign 
the uncertainties listed in the PDG [3TJ for the R — > 
transition helicity amplitudes. The main uncertainty is 
due to 5*11(1535), and we believe that this estimate of un- 
certainties is very conservative. The actual fit describes 
the data in the second resonance region certainly bet- 
ter than ±70%. We believe that this uncertainty can be 
further reduced. 

Finally, we obtain for the hadronic 0(a 5 ) contribution 
to the 2P — 2S Lamb shift in muonic deuterium set forth 
in Eq. ^ 



AE = (-40±5)/zeV. 



(36) 



V. DISCUSSION 

We have split the contribution of the nucleon's Comp- 
ton tensor to the Lamb shift of the muonic hydrogen 
atom into three parts, E el , E inel and E subt . The first 
two, corresponding to elastic scattering off the proton 
and photoexcitation of resonances are in agreement with 
previous work by other authors. The last term contains 
the contribution of the real subtraction to the Compton 
tensor and is the only one where significant uncertainty 



has remained. Specifically, in the analyses of [5] the sub- 
traction function was identified with 



T!(0,Q 2 ) = -^F 2 D (Q 2 ) + Q 2 (3(Q 2 ), 



(37) 



where Fd{Q 2 ) stands for the Dirac form factor, and 
(3(Q 2 ) for the generalized magnetic polarizability that for 
real photons reduce to the usual magnetic polarizability 
of Compton scattering, (3(0) = (3 m- Its Q 2 dependence 
was taken by analogy with elastic form factors. In Ref. 
[TT] it was argued that 



T\(0,Q 2 ) = - 



AI 



(38) 



where we put together the two contributions identified in 
[HI as Tf s (0,Q 2 ) = Q 2 (3(Q 2 ) and T^ no ~ pole (0,Q 2 ) = 
— jj for clarity. The common feature of the two approx- 
imations is that at Q 2 = they reduce to the Thomson 
term. However, they differ already in the first derivative, 
and they effectively operate with two different values of 
(3 that is a measured quantity. We define 



Ti(Q 2 



Ti(0,Q 2 ) + f 



Q 2 



(39) 



the function that enters the calculation of the Lamb shift, 
and evaluate this function at Q 2 = 0. With the model of 
2 one obtains 



Ti(0) 



M 



while the model of Ref. [TT] gives 
?i(0)=/3. 



(40) 



(41) 



The difference is not small and amounts to 3.3xl0~ 4 
fm 3 , of the same size as the polarizability itself. What 
complicates the issue is the impossibility to measure 
Ti(0, Q 2 ) directly since the kinematical arguments are 
in the unphysical region. The problem of low-energy ex- 
pansion of doubly virtual Compton scattering was ap- 
proached by two of us in [21] in terms of a fully model- 
independent low-energy theorem. It was found that it is 
only possible to unambiguously identify Ti(0,Q 2 ) with 
a combination of known or measurable quantities (form 
factors and polarizabilitics) modulo a dispersion inte- 
gral in the annihilation channel that is largely unknown. 
Rewriting the findings of Ref. [UJ for Ti(0, Q 2 ) we find 



Ti(0,Q 2 



AI 



[F 2 D {Q 2 ) 



tF 2 p (Q 2 )} 



Q 2 m 2 



(42) 



where we omitted terms coming from that dispersion in- 
tegral in the annihilation channel. The reason for such 
detailed discussion is to remind the reader that to relate 
the unphysical subtraction constant Ti(0, Q 2 ) to measur- 
able quantities like the polarizability and elastic form fac- 
tors, a good deal of caution should be exercised. 
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Carlson & Vanderhaeghen 



Pachucki 




Q z (GeV 2 ) 



FIG. 4: (Color online) Subtraction function [7\(0,<2 2 ) - 
Ti(0,0)]/Q 2 in units of 1(T 4 fm 3 as obtained from FESR 
(solid), from the model of Ref. [11] (dashed) and from Ref. 
[2] (dash-dotted). 



Following the analysis presented in this paper, the sys- 
tematic uncertainty in the Lamb shift from this term 
has been significantly reduced. We have employed the 
method of the Finite Energy Sum Rules to analyze this 
term, explicitly displaying the contributions it receives 
from the known i-channel Regge and s-channel reso- 
nances. There is no double counting of these resonances 
with respect to E mel . The alternative analysis presented 
here provides information on the subtraction term from 
Regge theory and the resonance region, reducing the un- 
knowns to the fixed pole of Compton scattering. Our 
Finite Energy Sum Rule in Eq. ( |30| has for the first 
time made it possible to predict the Q 2 -dependence of 
the subtraction function directly from existing experi- 
mental data. In Fig. [4] we compare the function Ti(Q 2 ) 



as obtained from FESR to phenomenological Ansatze of 
previous analyses. We observe that all approaches effec- 
tively have similar values of 7\ (0) but in view of the com- 
plicated situation with the low-energy theorem discussed 
above we stress that this is a coincidence. Neglecting 
the i-channel contributions in Eq. (42) and removing 



the contributions of the form factors (373 xlO -4 fm 3 and 
1.5xl0 -4 fm 3 ) we would arrive at /? = —1.8 x 10~ 4 fm 3 . 

We have shown that the contribution of the subtrac- 
tion term AE subt is small, « 3^eV, and its large relative 
error of order 5/ieV does not alter the conclusion that 
the overall contribution of the nucleon photoexcitation 
processes to the Lamb shift in muonic hydrogen is about 
-40±5/ieV. This is in agreement with the recent eval- 
uation of Carlson and Vanderhaeghen (— 37 ± 2.5/ieV), 
both being somewhat larger than earlier determinations 
of order — 20/ieV. Our overall estimated uncertainty has 
increased a bit respect to earlier work [5] as well as chi- 
ral perturbation theory [35], but we feel we have better 
control of systematic unknowns. 

The 300/z eV discrepancy between the direct muonic 
Hydrogen Lamb shift measurement and estimates for it 
based on usual (electronic) Hydrogen is unnaturally large 
for the hadronic structure-dependent corrections at order 
0(a 5 ) that have been proposed in the literature, basi- 
cally Eq. and the explanation must be looked for 
elsewhere. 
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